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Objectives Numerous studies have revealed the adverse health effects of acute and 
chronic exposure to particulate matter less than lo pm in aerodynamic diameter (PMio). 
The aim of the present study was to examine the spatial distribution of PMio concentra- 
tions and cardiovascular mortality and to investigate the spatial correlation between 
PMio and cardiovascular mortality using spatial scan statistic (SaTScan) and a regres- 
sion model. 

Methods From 2008 to 2010, the spatial distribution of PMio in the Seoul metropolitan 
area was examined via kriging. In addition, a group of cardiovascular mortality cases 
was analyzed using SaTScan-based cluster exploration. Geographically weighted regres- 
sion (GWR) was applied to investigate the correlation between PMio concentrations and 
cardiovascular mortality. 

Results An examination of the regional distribution of the cardiovascular mortality was 
higher in provincial districts (gu) belonging to Incheon and the northern part of Gyeonggi- 
do than in other regions. In a comparison of PMio concentrations and mortality cluster 
(MC) regions, all those belonging to MC 1 and MC 2 were found to belong to particulate 
matter (PM) 1 and PM 2 with high concentrations of air pollutants. In addition, the GWR 
showed that PMio has a statistically significant relation to cardiovascular mortality. 
Conclusions To investigate the relation between air pollution and health impact, spatial 
analyses can be utilized based on kriging, cluster exploration, and GWR for a more sys- 
tematic and quantitative analysis. It has been proven that cardiovascular mortality is 
spatially related to the concentration of PMio. 
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Introduction 

Numerous studies on environmental epidemiology have re- 
ported the adverse effects of short- and long-term exposure to 
air pollution. Such exposure can exacerbate pre-existing heart 
and lung diseases^ and it can lead to cardiovascular disease^ re- 
spiratory disease; asthma^ impaired lung function^ and prema- 
ture death [l;2]. 

Particulate matter less than 10 \\m in aerodynamic diameter 



(PMio) increases oxidative stress and inflammatory response^ 
and it is also associated with autonomic nervous system disor- 
ders and changes in blood coagulation^ which can lead to car- 
diovascular disease and cardiovascular-related mortality [3-5]. 

Environmental epidemiology have quantitatively evaluated 
the adverse health effects caused by air pollution through vari- 
ous statistical analyses^ including time-series studies and case- 
control studies [6;7]. Studies that incorporate the geospatial 
analysis method are actively being conducted. The geographic 
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information system can effectively analyze spatial information 
by managing and analyzing spatial data obtained through a pre- 
cise and systemic review of the physical space [8]. Spatial analy- 
sis has become an effective tool in the field of environmental ep- 
idemiology [9]. 

Health levels vary between different regions^ and thus in the 
field of environmental health; it is important to identify regional 
characteristics and distribution. Health effects due to exposure 
to environmental pollutants are subject to spatial analysis^ and it 
is necessary to determine the spatial correlation between them 
[10]. 

In a spatial epidemiology spatial analyses such as measuring 
distance; identifying clusters; spatial smoothing and interpola- 
tion; and spatial regression were used; and such methods are 
widely used when researching health effects with respect to air 
pollution [ll]. To identif)^ the spatial correlation between po- 
tential health effects and environmental hazards; it is useful to 
determine cluster of health effects; assess the spatial distribution 
of risk for disease; and conduct a spatial analysis comparing en- 
vironmental data with health data [ 12; 1 3] . 

In this study the spatial analyses were used to determine the 
correlation between PMio and cardiovascular mortality in the 
Seoul metropolitan area. To identify the characteristics of spatial 
distribution with regard to PMio and cardiovascular mortality 
spatial analysis was conducted in each location; and cluster anal- 
ysis was conducted to examine areas with high a cardiovascular 
mortality rate. Moreover, in order to determine the spatial cor- 
relation between PMio and cardiovascular mortality spatial dis- 
tribution of PMio and cluster areas of cardiovascular mortality 
were compared; and a spatial correlation was assessed quantita- 
tively using regression analysis. 

Materials and Methods 

Study Scope and Data 

The study areas were 79 provincial districts in Seoul; Incheon; 
and Gyeonggi-do within the metropolitan area (25 in Seoul; 1 1 
in Incheon; and 44 in Gyeonggi-do). The metropolitan area is 
11;745 km^; covering approximately 11.8% of South Koreas 
land area (99;720 km^). As of 2010; the population of Seoul 
was 9;794;000; that of Incheon was 2;663;000; and that of 
Gyeonggi-do was 113;790;000. These three areas account for 
49.1% of the entire South Korean population. 

As of 2010; the National Institute of Environmental Research 
has been operating 102 urban air monitoring stations in the 
metropolitan area in order to assess the level of air quality in the 
city. Data from the urban air monitoring network were used for 



PMio; and the ^-ray absorption method is being used to mea- 
sure PMio on an hourly basis. 

We obtained daily counts of deaths between January 1; 2008 
and December 31; 2010 from the National Statistical Office; 
South Korea. This study focused on deaths caused by cardiovas- 
cular-related diseases (international Classification of Diseases 
10th revision [ICD-IO]; code 100-199) within the metropolitan 
area; based on the time of death. 

Spatial Distribution of PMio 

To investigate the spatial distribution of PMio; kriging; which 
is widely used to estimate the spatial distribution of environ- 
mental factors; was applied. To estimate the concentration of 
PMio unobserved in this study ordinary kriging was applied by 
using the spherical theoretical variogram model as the weight. 
Compared to the Gaussian model or the exponential model; the 
calculated error value is minimal in the spherical model [14]. In 
general; the formula for ordinary kriging is as follows: 

i=i 

Here; z* is the critical value of the given point; Zi is the value 
from the established point; ji is the weight for each of the neigh- 
boring data used; and n is the number of data used for the krig- 
ing calculation. Weight is a function of distance; and it should be 
determined so that there is minimal variance of estimation be- 
tween the predicted value and the observed value. 

To assess the concentration of PMio in the Seoul metropolitan 
area between 2008 and 2010; the concentration of PMio in the 
79 provincial districts was calculated using kriging and zonal 
statistics via ArcGIS 10.1 (ESRI Inc.; Redlands; CA; USA). 

Software for The Spatial and Space-Time Scan Statistic 

To assess the overall spatial correlation; Morans index (l) was 
calculated for the distribution of cardiovascular mortality [15]. 
Morans I ranges from -1 to +1; with a positive /negative sign 
representing positive/negative spatial autocorrelation [16]. 
OpenGeoDa was used to determine the spatial autocorrelation 
of cardiovascular mortality. 

Cluster detection was performed using software for the spatial 
and space-time scan statistic (SaTScan) version 9.0 (Martin 
KuUdorff; Boston; MA; USA) with a Poisson model. Statistical 
cluster areas of hot spots are mapped based on SaTScan s likeli- 
hood ratio [17]. This study used SaTScans Poisson model; and 
the likelihood ratio y of the Poisson model is as follows: 
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Here; C is the entire value and c is the observed value. E[c] is 
the predicted value of the observed value and J() is the index 
factor. 

Geographically Weighted Regression 

Geographically weighted regression (GWR) was performed in 
order to analyze the spatial correlation between PMio and car- 
diovascular mortality Compared to general regression analysis^ 
which represents the study area with one regression coefficient^ 
GWR estimates regression coefficients regionally thereby ex- 
plaining the changes in the independent variable depending on 
the location of the spatial unit [ 18] . 

GWR uses the weighted matrix (W) between regions to calcu- 
late l^jif each of the regions regression coefficients^ with regard to 
variable;'; which depends on the geographical location /. GWRs 
regression analysis is as follows: 

) = 1 

Jji=(X'WiXy\X'WiY) 

Here; / is the location of an area; Y represents cardiovascular 
mortality the independent variable is the PMio concentration; 
and the regression coefficients and error term are represented as 

and e, respectively. The geographical weighted estimate of re- 
gression coefficient (jS/j) is calculated by the weighted least 
square; and W, is the local degree of the neighborhood / [19]. 
The GWR used ArcGIS version 10.1 to determine the relation- 
ship between PMio and cardiovascular mortality 

Results 

Table 1 displays the basic statistics of the number of cardiovas- 
cular mortality and PMio concentration from 2008 to 2010; 
which were identified through kriging. In 2010; the annual aver- 
age concentration of PMio was 52.06 [xg/m^; showing consider- 
able improvement in contrast to 2008 and 2009. The number of 



cardiovascular mortality was maintained at a constant level dur- 
ing the research period; and the number of cases over three 
years was 778.34. 

The global Morans I showed that the reported rates were 0.36 
in 2008; 0.29 in 2009; 0.46 in 2010; and 0.44 between 2008 and 
2010; all of which proved to be statistically significant. The Mo- 
rans I was 0.2 or higher; indicating that the cardiovascular mor- 
tality presented a spatial autocorrelation; and in particular; Mo- 
rans I was higher than 0.4 in the case of cardiovascular mortality 
in 2010; along with the period between 2008 and 2010; signify- 
ing a strong spatial autocorrelation. 

Standardized mortality ratio (SMR) and the results of spatial 
scan statistics for cardiovascular mortality are shown in Figure 1. 
The SMR for each provincial district (gu) of the Seoul metro- 
politan area was 0.55-1.38 in 2008; 0.51-1.48 in 2009; and 
0.70-1.69 in 2010. The SMR was higher within Incheon and 
the northern part of Gyeonggi-do compared to other regions 
between 2008 and 2010. In contrast; the SMR of Seoul and the 
neighboring Gyeonggi-do area was lower than that in other re- 
gions. 

The difference in the likelihood ratio determined the cluster 
region; and mortality cluster (MC) 1 indicates the likelihood ra- 
tio within the highest risk cluster. Regions MC 1 and MC 2 
were categorized as the provincial district within Incheon and 
the northern part of Gyeonggi-dO; respectively. 

Table 2 displays the results of SaTScan such as the population 
exposed to risk within the cluster as a hot spot for cardiovascular 
mortality between 2008 and 2010; incidences of cardiovascular 
mortality as well as relative risk; and data regarding the cluster, 
which were statistically significant. 

The concentration of PMio was categorized in four equal in- 
tervals in order to compare the distribution of regional clusters 
regarding PMio and cardiovascular mortality. Particulate matter 
(PM) 1 signifies a region with PMio concentration that is higher 
than 75%; PM 2 50-75%; PM 3 25-50%; and PM 4 lower than 
25%. Figure 2 indicates overlapping spatial distribution regions 



Table 1. Summaries of PMio and cardiovascular death in the Seoul metropolitan area of South Korea, 2008-2010 







Minimum 


25th 


Percentiles 
50th 


75th 


Maximum 


Mean 


Standard 
deviation 


PMio (|ig/m^) 


2008 


49.32 


54.09 


56.63 


58.52 


61.54 


56.29 


2.67 




2009 


51.28 


54.14 


56.53 


59.84 


63.76 


57.21 


3.49 




2010 


43.55 


48.99 


52.25 


55.91 


59.03 


52.06 


4.32 




2008-2010 


48.41 


52.36 


55.22 


57.94 


61.12 


55.19 


3.29 


Cardiovascular death 


2008 


26 


168 


241 


359 


564 


262.47 


118.48 




2009 


34 


162 


236 


348 


526 


252.77 


112.04 




2010 


42 


182 


240 


339 


539 


263.10 


111.31 




2008-2010 


34 


180 


235 


345 


536 


259.45 


112.58 



PMio, particulate matter less than 10 pm in aerodynamic diameter. 
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Figure 1. Standardized mortality ratio (SMR) (A) and mapping of cardiovascular mortality clusters (MC) (B) with higher mortality rates from cardiovascular 
disease. RR, relative risk. 



Table 2. Significant clusters found by the spatial scan statistic software for cardiovascular mortality in the Seoul metropolitan area of South Korea, 2008-201 0 



Year 


Cluster label (n) 


Population at risk 


No. of case 


Relative risk 


Log likelihood ratio 


p-value 


2008 


MC 1 (22) 


7,456,536 


6,774 


1.17 


55.68 


<0.001 




MC 2 (3) 


268,308 


472 


1.25 


11.06 


0.001 


2009 


MC 1 (8) 


2,561,111 


2,057 


1.26 


56.73 


<0.001 




MC 2 (3) 


277,285 


521 


1.38 


23.96 


<0.001 




MC 3 (2) 


477,862 


484 


1.25 


10.78 


0.001 


2010 


MC 1 (8) 


2,561,111 


2,613 


1.27 


60.25 


<0.001 




MC 2 (4) 


318,580 


600 


1.26 


14.65 


<0.001 




MC 3 (5) 


1,127,049 


1,265 


1.16 


12.73 


<0.001 


2008-2010 


MC 1 (8) 


2,561,111 


7,683 


1.26 


166.86 


<0.001 




MC 2 (4) 


318,580 


1,780 


1.27 


46.18 


<0.001 




MC 3 (5) 


1,127,049 


3,660 


1.14 


27.34 


<0.001 



MC, mortality cluster. 



that are the concentration range of PMio and hot spots of car- 
diovascular mortality. 

Table 3 presents the relationships between MCs and PMio con- 
centration. The results of a comparison of the average PMio 
range and the cluster as hot spots of cardiovascular mortality ac- 
cumulated throughout 2008-2010 show that the eight provincial 
districts of the MC 1 area all coincided with the PM 1 area^ 
which are regions with a high concentration of PMio. A compari- 
son of the corresponding ratio regarding the cardiovascular MC 



regions by the level of PMio concentration indicates that cluster 
of cardiovascular mortality correspond mostly with areas of high 
PMio concentration. 

In order to quantitatively analyze the correlation between the 
SMR of cardiovascular mortality and PMio; spatial regression 
analysis using the GWR was conducted. This study excluded 
the 1% region with the highest level of PMio concentration (one 
provincial district) as an outlier among the 79 provincial dis- 
tricts of the Seoul metropolitan area. 
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I 1 43.56-52.99 53.00-55.33 ^^55.34-58.18 




2008 



2009 



2010 



2008-2010 



Figure 2. Spatial distribution of PIVlio concentrations and cardiovascular mortality clusters (MC). PMio, particulate matter less than 10 pm in aerodynamic 
diameter; RR, relative risk. 



Table 3. Relationship between cardiovascular mortality clusters (MC) and 
PMio concentration 



MC 



PMio concentration 



2008 


PM 1 (26) 


PM2(21) 


PM 3 (22) 


PM4(10) 


MC 1 (22) 


12(54.5) 


10(45.5) 






MC 2 (3) 


1 (33.3) 


2 (66.7) 






2009 


PM 1 (29) 


PM2(18) 


PM 3 (28) 


PM 4 (4) 


MC 1 (8) 


8(100.0) 








MC 2 (3) 


1 (33.3) 


2 (66.7) 






MC 3 (2) 






1 (50.0) 


1 (50.0) 


2010 


PM 1 (4) 


PM2(21) 


PM 3 (9) 


PM 4 (45) 


MC 1 (8) 


3 (37.5) 


5 (62.5) 






MC 2 (4) 




4(100.0) 






MC 3 (5) 




3 (60.0) 


2 (40.0) 




2008-2010 


PM1 (17) 


PM 2 (22) 


PM3(16) 


PM 4 (24) 


MC 1 (8) 


8(100.0) 








MC 2 (4) 


1 (25.0) 


3 (75.0) 






MC 3 (5) 


1 (20.0) 


4 (80.0) 







Values are presentd as number (%). 

PMio, particulate matter less than 10 pm in aerodynamic diameter; PM, parti- 
culate matter. 



Table 4 presents the results of the correlation of cardiovascular 
mortality and PMio using the GWR; which signifies a significant 
relevance between PMio and cardiovascular mortality. The min- 
imum; maximum; and average regression coefficient of PMio 
and cardiovascular mortality by region in the GWR were 0.799; 
1.562; and 0.956; respectively and a difference in the correlation 
between PMio and cardiovascular mortality was reported. The 
R^fld; of the GWR was 0.205; and the result of analyzing the re- 
sidual of Morans I was 0.039; which indicates that the spatial 
autocorrelation of the residual shows an arbitrary pattern and 
satisfies the assumption of the basic GWR. 

Discussion 

In this research; various spatial epidemiological analyses were 



Table 4. Geographically weighted regression (GWR) of PMio and cardiovas- 
cular mortality in the Seoul metropolitan area of South Korea, 2008-2010 

GWR 



Min p Max p Mean p 



SE 



PMio 
Local R2 

R\dj 

AlC 

Moran's index of residual 



0.799 
0.157 



1.562 0.956 
0.249 0.198 

0.205 
639.324 

0.039** 



0.102 
0.017 



Min, minimum; Max, maximum; SE, standard error; PMio, particulate matter 
less than 10 pm in aerodynamic diameter; AlC, Akaike information criterion. 
**p<0.01 . 



applied to the Seoul metropolitan area between 2008 and 2010 
in order to determine the correlation between PMio and cardio- 
vascular mortality. Kriging was performed on PMio concentra- 
tion for to determine spatial distribution; and spatial autocorre- 
lation was determined by using Morans I. Disease mapping of 
cardiovascular mortality was conducted with use of SaTScan in 
order to explore the geographical distribution; and cluster of 
cardiovascular mortality were investigated. AlsO; the spatial cor- 
relation between PMio and cardiovascular mortality was identi- 
fied quantitatively using the GWR. 

Disease surveillance researches pursuing the temporal spatial 
distribution of disease using SaTScan have been actively con- 
ducted. Fukuda et al. [20] have investigated clusters of mortality 
cases due to colon cancer; lung cancer; and breast cancer using 
SaTScan throughout Japan as the subject matter; and the rela- 
tion between mortality clusters and social clusters (socio-eco- 
nomic indicators and population density) was analyzed. The re- 
sults of the analysis revealed that the mortality cluster due to 
lung cancer and breast cancer deaths were the urban areas. 
Horst and Coco [21] used SaTScan to analyze clusters of pa- 
tients visiting basic health services and emergency rooms due to 
respiratory and gastrointestinal diseases. Su et al. [22] used 
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SaTScan to spatially analyze clusters of cancer patients based on 
gender; type of cancer^ and ethnicity in order to identify the 
range and characteristics of cancer center service areas. 

AlsO; the following five categories were examined for the pur- 
pose of disease surveillance: data processings analysis methods^ 
technical issues^ analysis output^ and user s facility regarding the 
four spatio temp oral analysis programs (SaTScan^ ClusterSeer^ 
GeoSurveillance^ and R- Surveillance); results showed that 
SaTScan was the most robust program for searching clusters of 
health effects [23]. 

Disease surveillance entails collecting; analyzing; and inter- 
preting data in order to identif)^ the frequency and distribution 
of diseases within the population [24]. It is an important pro- 
cess that supplies valuable information for policy decisions and 
implementation. Utilizing spatial cluster analysis enables basic 
data to be used for preventative activities by identifying areas 
with relatively high disease risk as well as disease clusters. 

Internationally research analysis on the relationship between 
air pollution and health effects is being conducted using various 
spatial analyses; and a significant spatial correlation between 
health and air pollution is being reported. Corburn et al. [25] 
used SaTScan to identif)^ asthma cluster in New York City and 
they analyzed the correlation between land use; housing type; 
and air pollution. The analysis results revealed a quantitative 
correlation of 69% between asthma cluster and air pollution fa- 
cilities. Jephcote and Chen [26] utilized the GWR in order to 
examine the relationship between respiratory hospitalizations of 
children under the age of 15; socio-economic indicators; and 
PMio from automobile emissions in Leicester, UK. The results 
of their study indicated that respiratory hospitalization of chil- 
dren is related to socio-economic indicators and PMio from au- 
tomobile emissions; and this relationship differed between re- 
gions. In this research; a comparison of the PMio and the hot 
spot of cardiovascular mortality identified through SaTScan 
showed that in most caseS; the cardiovascular mortality cluster 
was included within the area with high levels of PMio. AlsO; a 
significant correlation was identified as a result of analyzing the 
PMio concentration and cardiovascular mortality using the 
GWR; and the spatial correlation appeared differently in the re- 
sult of the regional regression coefficient calculation; which in- 
dicates a similar result as the previous research results. The rea- 
son that the regional spatial correlation between PMio and car- 
diovascular mortality appears differently is due to various fac- 
tors; such as concentration levels and components of the PMio; 
demographic characteristics; and sensitivity of the population 
exposed to the PMio; as well as socio-economic factors [27;28]. 

The spatial analysis using the GWR scientifically suggests re- 
gional differences between air pollution and health effects; and 



conducting policy intervention based on this information can 
help in the effective distribution of limited resources by priori- 
tizing the policy statement and the target area. 

Through the development of a geographic information system 
and the joining of statistical techniques; a more scientific analy- 
sis regarding the relationship between health risks due to envi- 
ronmental exposure and risks within the spatial epidemiology is 
made possible. Spatial epidemiology involves integrated analy- 
sis; such as visualization; exploratory analysis; and modeling 
based on spatial data regarding environmental exposure and 
health effects; and it examines the analysis of the spatial correla- 
tion and patterns in order to predict and prevent health risk 
[10;29]. This research utilized distance interpolation; disease 
mapping; and GWR analysis among various spatial methods in 
order to identify the correlation between air pollution and 
health effects; and it produced significant results of the correla- 
tion between air pollution and health effects through systematic 
and quantitative analysis. 

However; there are certain limitations to this research. First; it 
is an ecological study that uses population data; and it does not 
integrate personalized data due to the limitations of the data. 
Therefore; the research does not consider personal characteris- 
tics or dynamic factors such as population movement because 
the analysis is based on population data. Second; cardiovascular 
mortality can occur due to factors other than air pollution; in- 
cluding environmental factors; socio-economic factors; medical 
factors; and dietary factors. This study only analyzed health ef- 
fects caused by PMio. In the future; it will be necessary to con- 
figure the model to consider a number of factors. 

The spatial analyses and findings of the research on the rela- 
tionship between air pollution and health effects can be utilized 
as a scientific method to identif)^ vulnerable areas for evidence- 
based public health policies. 
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